Viriato: a Fourier-Hermite spectral code for strongly magnetised fluid-kinetic plasma dynamics

We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] (which reduce to the standard Reduced-MHD equations in the appropriate limit) and (ii) the kinetic reduced MHD (KRMHD) equations [Schekochihin et al., Astrophys. J. Suppl. 182:310 (2009)]. Two main applications of these equations are magnetised (Alfvenic) plasma turbulence and magnetic reconnection. Viriato uses operator splitting (Strang or Godunov) to separate the dynamics parallel and perpendicular to the ambient magnetic field (assumed strong). Along the magnetic field, Viriato allows for either a second-order accurate MacCormack method or, for higher accuracy, a spectral-like scheme composed of the combination of a total variation diminishing (TVD) third order Runge-Kutta method for the time derivative with a 7th order upwind scheme for the fluxes. Perpendicular to the field Viriato is pseudo-spectral, and the time integration is performed by means of an iterative predictor-corrector scheme. In addition, a distinctive feature of Viriato is its spectral representation of the parallel velocity-space dependence, achieved by means of a Hermite representation of the perturbed distribution function. A series of linear and nonlinear benchmarks and tests are presented, including a detailed analysis of 2D and 3D Orszag-Tang-type decaying turbulence, both in fluid and kinetic regimes.Comment: 42 pages, 15 figures, submitted to J. Comp. Phy

Comprehensive evaluation of the linear stability of Alfvén eigenmodes driven by alpha particles in an ITER baseline scenario

The linear stability of Alfvén eigenmodes in the presence of fusion-born alpha particles is thoroughly assessed for two variants of an ITER baseline scenario, which differ significantly in their core and pedestal temperatures. A systematic approach based on CASTOR-K (Borba and Kerner 1999 J. Comput. Phys. 153 101; Nabais et al 2015 Plasma Sci. Technol. 17 89) is used that considers all possible eigenmodes for a given magnetic equilibrium and determines their growth rates due to alpha-particle drive and Landau damping on fuel ions, helium ashes and electrons. It is found that the fastest growing instabilities in the aforementioned ITER scenario are core-localized, low-shear toroidal Alfvén eigenmodes. The largest growth-rates occur in the scenario variant with higher core temperatures, which has the highest alpha-particle density and density gradient, for eigenmodes with toroidal mode numbers . Although these eigenmodes suffer significant radiative damping, which is also evaluated, their growth rates remain larger than those of the most unstable eigenmodes found in the variant of the ITER baseline scenario with lower core temperatures, which have and are not affected by radiative damping

Sensitivity of alpha-particle-driven Alfvén eigenmodes to q-profile variation in ITER scenarios

A perturbative hybrid ideal-MHD/drift-kinetic approach to assess the stability of alpha-particle-driven Alfv�n eigenmodes in burning plasmas is used to show that certain foreseen ITER scenarios, namely the MA baseline scenario with very low and broad core magnetic shear, are sensitive to small changes in the background magnetic equilibrium. Slight variations (of the order of ) of the safety-factor value on axis are seen to cause large changes in the growth rate, toroidal mode number, and radial location of the most unstable eigenmodes found. The observed sensitivity is shown to proceed from the very low magnetic shear values attained throughout the plasma core, raising issues about reliable predictions of alpha-particle transport in burning plasmas

Spatial resolution on a small animal RPC-PET prototype operating under magnetic field

It has been demonstrated in previous work that the RPC-PET technology is able to deliver radioisotope image resolutions approaching the physical limits of the PET principle. Here we study, by simulation, the effect of the magnetic field on the positron range to evaluate whether the spatial resolution of the RPC-PET could be improved by using an intense magnetic field. Six positron emitters of interest to small animal PET imaging (18F, 11C, 15O, 68G, 62Cu and 86Y) are considered. Results suggest that a three-fold improvement on the spatial resolution may be obtained under a magnetic field of 10 T for the higher energy radioisotopes like 86Y or 62Cu, and by about 20% for the lower energy ones, like 18F or 11C.http://www.sciencedirect.com/science/article/B6TVD-4KV4MR8-15/1/efa4356ba16a85fd7740937bd0e672a

Very high position resolution gamma imaging with resistive plate chambers

In this study we present experimental results from a first prototype of a positron emission tomography system based on the resistive plate chamber (RPC) technology. The system is composed of two counting heads, each one containing 16 single-gap RPC detectors capable of detecting the photon interaction point in the transaxial plane. Uniformity studies were performed for image resolution and sensitivity, yielding a rather uniform image resolution close to 0.3 mm FWHM across the field of view. The contribution of the photons noncolinearity effect to the intrinsic spatial resolution was also studied, causing a variation from 0.52 to 0.63 mm when the system diameter ranges from 60 to 120 mm, in agreement with calculations.http://www.sciencedirect.com/science/article/B6TJM-4K5HV7F-9/1/7b655e08199d2b28acc8050106c2271

Performance of distributed multiscale simulations

Multiscale simulations model phenomena across natural scales using monolithic or component-based code, running on local or distributed resources. In this work, we investigate the performance of distributed multiscale computing of component-based models, guided by six multiscale applications with different characteristics and from several disciplines. Three modes of distributed multiscale computing are identified: supplementing local dependencies with large-scale resources, load distribution over multiple resources, and load balancing of small- and large-scale resources. We find that the first mode has the apparent benefit of increasing simulation speed, and the second mode can increase simulation speed if local resources are limited. Depending on resource reservation and model coupling topology, the third mode may result in a reduction of resource consumption

Performance of distributed multiscale simulations

Multiscale simulations model phenomena across natural scales using monolithic or component-based code, running on local or distributed resources. In this work, we investigate the performance of distributed multiscale computing of component-based models, guided by six multiscale applications with different characteristics and from several disciplines. Three modes of distributed multiscale computing are identified: supplementing local dependencies with large-scale resources, load distribution over multiple resources, and load balancing of small- and large-scale resources. We find that the first mode has the apparent benefit of increasing simulation speed, and the second mode can increase simulation speed if local resources are limited. Depending on resource reservation and model coupling topology, the third mode may result in a reduction of resource consumption. © 2014 The Authors

Modelling of the effect of ELMs on fuel retention at the bulk W divertor of JET

Effect of ELMs on fuel retention at the bulk W target of JET ITER-Like Wall was studied with multi-scale calculations. Plasma input parameters were taken from ELMy H-mode plasma experiment. The energetic intra-ELM fuel particles get implanted and create near-surface defects up to depths of few tens of nm, which act as the main fuel trapping sites during ELMs. Clustering of implantation-induced vacancies were found to take place. The incoming flux of inter-ELM plasma particles increases the different filling levels of trapped fuel in defects. The temperature increase of the W target during the pulse increases the fuel detrapping rate. The inter-ELM fuel particle flux refills the partially emptied trapping sites and fills new sites. This leads to a competing effect on the retention and release rates of the implanted particles. At high temperatures the main retention appeared in larger vacancy clusters due to increased clustering rate